Rate limit for photoassociation of a Bose-Einstein condensate 
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We simulate numerically the photodissociation of molecules into noncondensate atom pairs that 
accompanies photoassociation of an atomic Bose-Einstein condensate into a molecular condensate. 
Such rogue photodissociation sets a limit on the achievable rate of photoassociation. Given the atom 
density p and mass m, the limit is approximately 6hp 2 ^ 3 /m. At low temperatures this is a more 
stringent restriction than the unitary limit of scattering theory. 
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Photoassociation of a Bose-Einstein condensate 
has been predicted to feature examples of quantum co- 
herence, such as Rabi oscillations between atomic and 
molecular condensates ]||| , rapid B and Raman pj adi- 
abatic passage from an atomic condensate to a molecu- 
lar condensate, and a superposition of macroscopic num- 
bers of atoms and molecules || . The Feshbach resonance 
in a magnetic field is mathematically equivalent to pho- 
toassociation, and should give rise to analogous phenom- 
ena H . Spontaneous-emission decay of the photoassoci- 
ated molecules causes experimental problems JTJ , and so 
far experiments on photoassociation of a Bose-Einstein 
condensate have not demonstrated coherent optical 
transients. Nonetheless, the prospect of a molecular con- 
densate continues to motivate further work. 

Many of the past theoretical studies of coherent pho- 
toassociation have taken into account only the atomic 
and molecular condensates. Notably ignored is the pho- 
todissociation of a molecule into a pair of atoms, nei- 
ther of which belongs to the atomic condensate. We 
have investigated such "rogue" photodissociation [^]]7]] 
using essentially dimensional arguments, and suggested 
that there is a shortest possible time scale and a cor- 
responding fastest possible Rabi frequency for coherent 
photoassociation. Given the laser intensity /, the Rabi 
frequency fl is proportionally to the laser field strength, 
fl cx \fl, whereas the rate of rogue photodissociation 
scales as r cx /. The equality = T occurs when both 
of these frequencies are approximately equal to 
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Here p = pa + ^Pm is an invariant density derived from 
the densities of the atoms and the molecules, and p is 
the reduced mass of a pair of atoms. In an attempt to 
speed up coherent photoassociation beyond this limit by 
increasing the laser intensity, damping due to rogue pho- 
todissociation overtakes coherent photoassociation 0,]?). 
There have also been numerical case studies on the effects 
of noncondensate modes on photoassociation , but 

these do not spell out time or frequency scales compara- 
ble to Eq. (TO). 



In the present paper we formulate a minimal yet real- 
istic, microscopic model for rogue photodissociation, and 
solve it numerically. The results demonstrate that there 
is a maximum achievable photoassociation rate, of the 
order 6uj p . The existence of the rate limit is at variance 
with both the wave picture of condensate atoms, and the 
standard scattering theory. 

We take a plane-wave laser field with photon mo- 
mentum q to drive photoassociation and photodissoci- 
ation. Initially there is only a condensate of N zero- 
momentum atoms present. By momentum conservation, 
only molecules with momentum q will be generated in 
the primary photoassociation. In photodissociation such 
molecules break up into a pair of atoms with opposite 
momenta. Processes in which both atoms return to the 
condensate are the subject of our earlier work on coher- 
ent photoassociation [gj^], whereas here the focus is on 
processes in which the molecules break up into a pair of 
noncondensate atoms. 

Our description starts from our generic Hamiltonian 
for photoassociation and dissociation HQ . We write the 
Heisenberg equations of motion for the boson operators 
Qk (&k) characterizing atoms (molecules) with momen- 
tum k, and implement two main approximations. First, 
in the classical approximation analogous to the use of 
the Gross-Pitaevskii equation, we treat the atomic and 
molecular condensate operators a = ao and b = fe q as 
c-numbers. Second, we take into account all pairs of 
noncondensate atoms with opposite momenta, but no 
further photoassociation of noncondensate atoms to any 
molecular states except back to the molecular conden- 
sate. It then turns out that the equations of motion for 
the c- number atomic amplitude a — a/yfN, molecular 
amplitude /3 = y^2/N b, and noncondensate atom pair 
amplitudes Ck = (aka-k) close. As usual, we discuss s- 
wave photoassociation and dissociation, so that we may 
uniquely label the two-atom modes with their energy he. 
The equations of motion may be cast in the form 
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To explain the notation, we first consider the model in 
the absence of the coupling between condensate atoms 
and molecules oc fl. The quantity S is the detuning of 
the laser above the photodissociation threshold of the 
molecules (corrected for photon recoil energy and 
£(e) is the coupling strength for photodissociation of a 
molecule to an atom pair with energy he. Using the 
Fermi golden rule, we find the photodissociation rate for 
molecules T(S) — 27r|£((5)| 2 . In our modeling we turn the 
tables, and assert the equality 
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relating the matrix element £(e) and the on-shell pho- 
todissociation rate T(e) to an atom pair with energy he. 
Second, we know from our theory of coherent photoasso- 
ciation [^,0 that the Rabi frequency fl is related to the 
rate of photodissociation by 
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Third, in order for Eq. (Q) to make sense, for low en- 
ergies r(e) oc y/e must hold true. This, in fact, is 
known as the Wigner threshold law. In our modeling 
we take the threshold law to hold for all relevant ener- 
gies. Finally, for notational convenience, we write the 
continuum amplitude in terms of a dimensionless quan- 
tity C{x) = c X Qyi~ x / 2 . Our final equations of motion for 
the amplitudes read 
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8 = iSB + i— = or 
V2 



The parameter 



r(fi) 

2nfl 



1 

8^ 



3/2 



(5a) 

(5b) 
(5c) 

(6) 



characterizes the relative importance of rogue photodis- 
sociation and coherent photoassociation. 

In Eqs. (g) we could readily absorb the frequency fl 
into the unit of time, but we will not do so as this would 
cause contorted discussions below. Coupling to the non- 
condensate modes introduces the frequency to p . 

We study Eqs. (^|) numerically, by discretizing the dis- 
sociation continuum. We also cut off the continuum at 
a maximum frequency em- Basically, the vector made of 
the amplitudes a, 8 and C(x n ), %b, satisfies a Schrodinger 
equation of the form itp = H(ip)i/j. As Eqs. (^|) are non- 
linear, the Hamiltonian H depends on the vector tp itself. 



However, if this dependence is ignored, the resulting lin- 
ear Schrodinger equation may be solved using the norm 
conserving algorithm described in Ref. Moreover, 
the matrix of the Hamiltonian is such that the run time 
scales linearly with the number of the continuum states 
n c . We take care of the nonlinearity with a predictor- 
corrector method. To execute a time step t — ► t + h, 
we first integrate the Schrodinger equation for the con- 
stant Hamiltonian H[ip(t)] to obtain the "predicted" 
wave function i/)(t + h), and then "correct" by integrat- 
ing the Schrodinger equation anew using the Hamiltonian 
\{H[ij}(t)\ + H[ip(t + h)]}. The algorithm is norm con- 
serving, the error for integrating over a fixed time inter- 
val scales with the step size at least as 0(h 2 ) (empirically 
often better), and the run time is linear in the number 
of continuum states. Typically we discretize the contin- 
uum in n c — 10 4 steps with the maximum frequency 
ejvf = 100 fl, but when needed for convergence we do not 
hesitate to resort to, say, n c = 10 6 and cm = 10 3 fl. 

Coupling to the dissociation continuum also shifts the 
molecular level JlJJlIt ■ We correct for the shift by adding 
to the detuning the Stark shift calculated using perturba- 
tion theory for on-threshold photodissociation, namely, 
r(ejvf)/7T- The shift may be very large; in fact, if we did 
not truncate the continuum, it would be infinite. This is 
our version of the renormalization problem discussed in 
Ref. H). In the limit of a large dissociation rate, per- 
turbation theory becomes unreliable. Correspondingly, 
it gets difficult to tell exactly what the detuning is with 
respect to the photodissociation threshold. This prob- 
lem is not unique to our model, but would emerge in real 
experiments as well. 

For demonstration, Fig. 1 gives the quantities Pa = 
\a\ and Pm = \@\ 2 , the probabilities that an arbitrary 
atom in the system is part of the atomic or the molec- 
ular condensate, as a function of time, given that the 
system started as an atomic condensate. In Fig. 1(a) 
we set 7 = 0, and thereby disregard rogue photodisso- 
ciation. For this (renormalized) detuning 5 — fl, the 
conversion between atoms (dotted line) and molecules 
(dashed line) oscillates in such a way that half of the 
atomic condensate is periodically converted into a molec- 
ular condensate. In Fig. 1(b) we add rogue photodisso- 
ciation, setting r(fi) = fl and hence 7 = l/(27r). After 
a few damped oscillations, the probabilities settle down. 
The total probability that an atom is part of either con- 
densate (solid line) drops to 0.8, indicating that in the 
process 20 % of the atoms are lost to rogue photodissoci- 
ation. Although the couplings between the condensates 
and from the molecular condensate to the noncondensate 
atoms remain, the atomic and the molecular condensates 
lock up in a superposition immune to rogue photodisso- 
ciation. It is known from the past work on decay of a 
discrete state to a shaped continuum (e.g., photodisso- 
ciation of a negative ion that, given a threshold to 
the continuum, some population may stay permanently 
trapped in the bound state. Curiously, the same applies 
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here to a superposition of atomic and molecular ampli- 
tudes, even though the condensate system is nonlinear 
and the superposition principle does not hold. 

Our main subject is limitations on conversion from 
atoms to molecules due to rogue photodissociation. In 
the traditional perturbation theory one would study the 
photoassociation rate as a function of the problem pa- 
rameters. However, as is obvious from Fig. 1, atom- 
molecule conversion is not governed by a simple expo- 
nential time dependence. A clear-cut photoassociation 
rate generally does not exist. Because of level shifts, 
it may also be difficult to know the true value of the 
detuning. To circumvent these problems, we have de- 
veloped a protocol that mimics conceivable experiments. 
We first fix an interaction time, call it t, once and for 
all. Our basic numerical operation is to calculate the 
loss of condensate atoms as a function of the Rabi fre- 
quency f2, rogue photodissociation amplitude 7, and de- 
tuning 6,1 — -Pa(^, 1, S; t), starting from a pure atomic 
condensate. We vary the detuning S, and find what 
we interpret as the maximal photoassociation probabil- 
ity with the other parameters fixed, 1 — P^(H,j;t). 
The optimal rate for photoassociation is then defined as 

R(n,r,T) = [i-p A {n,r,T)]/T. 

To get a quantitative handle on the variation of the 
effective rate R we divide it by u> p , and plot in Fig. 2 the 
result as a function of the variable l/(ru> p ) for four differ- 
ent Rabi frequencies, £1 = 1/r (dotted line), 2/r (dash- 
dotted line), 4/r (dashed line), and 8/r (solid line). It 
appears that there is a maximum rate for photoassocia- 
tion of the order of 6uj p . 

We emphasize once more that our maximal photoasso- 
ciation rate measures how many atoms one can remove 
in a fixed interaction time. There really does not have 
to be a rate at all. We demonstrate this by plotting in 
Fig. H the time dependences of the molecular probability 
Pm (dotted line), the probability for condensate atoms 
Pa (dashed line), and the total condensate probability 
Pa + Pm (solid line) as a function of time. The fixed 
parameters ft — 4/r, 6 — — 9.27 17, and 7 = 1.27 are 
chosen in such a way that this is the time evolution used 
to obtain the maximum of the curve labeled Q = 4/r in 
Fig. ||. The laser is tuned well below the photoassocia- 
tion threshold. There is very little molecular probability, 
so that the molecules are a virtual intermediate state. 
Instead, we have a slowly damped two-photon flopping 
between the atomic condensate and noncondensate atom 
pairs. It is doubtful that such flopping survives the vari- 
ous sources of incoherence left out from our model, such 
as photoassociation of noncondensate atoms to noncon- 
densate molecules or collisions of noncondensate atoms. 
But we surmise that added loss mechanisms further hin- 
der photoassociation, and our upper limit for the rate 
6w p would probably hold in a more complete model. 

So far we have discussed the case when the molecu- 
lar state is not damped by spontaneous emission. Our 
modeling is most appropriate for two-color photoassocia- 
tion in the Raman configuration. However, we have also 



carried out realistic simulations of one-color photoasso- 
ciation, in which the spontaneous decay rate of the pho- 
toassociated molecules is the largest frequency param- 
eter in the problem. Photoassociation once more does 
not proceed exponentially and there is some ambiguity 
in the definition of the rate, but the conclusion is that 
the maximum achievable rate for photoassociation very 
nearly equals ui p Jll| . 

If one is to think of atoms as particles, there must be 
a limit to photoassociation rate. In order to photoas- 
sociate, two atoms have to drift close to one another. 
This process takes the longer, the more dilute the gas is; 
and it cannot be sped up arbitrarily by increasing the 
light intensity. Beyond a certain intensity, the rate will 
no longer increase with increasing intensity. In fact, for 
a quantum mechanical, noninteracting, zero-temperature 
condensate, one can construct (up to a numerical factor) 
only one characteristic frequency, u) p . Not surprisingly, 
(up to a numerical factor) the maximum photoassocia- 
tion rate equals u> p 

On the other hand, within the two-mode approxima- 
tion, a gas can be photoassociated arbitrarily fast; just 
set 7 = in Eqs. (|5|) and let £1 — > 00 by increasing the 
laser intensity. In the two-mode approximation the atoms 
are waves, simultaneously everywhere in the atomic con- 
densate. The existence of a maximum to the photoas- 
sociation rate thus is a question of particle versus wave 
nature of the atoms. 

Interestingly, in scattering theory the argument goes 
in the opposite direction. Because of the wave nature of 
particles, the scattering cross section may be much larger 
than the geometric size of the colliding objects. The max- 
imum quantum mechanical cross section for photoassoci- 
ation, the unitary limit, gives the same rate limit as rogue 
photodissociation when the temperature is of the order of 
the critical temperature for Bose-Einstein condensation. 
However, at lower temperatures the newly found particle 
nature dominates, and the rate limit due to rogue pho- 
todissociation is more stringent than the unitary limit. 
There is no rogue-photodissociation limit in the usual res- 
onance scattering analysis of photoassociation |l7| , and 
an experimental observation would constitute a dramatic 
failure of scattering theory. 

The two-atom position correlation function straddles 
particle and wave pictures. Photoassociation evidently 
depletes two-body correlations at short distances, and 
the motion of the atoms has to replete them. Though 
we have not been able to make an explicit connection, 
we surmise that Ref. ]Tl| ] focusing on atom-atom corre- 
lations boils down essentially to the same physics as our 
rogue photodissociation. In Ref. |l(]] the Gross-Pitaevskii 
equations were actually integrated numerically in the mo- 
mentum representation. The phenomenon thus found 
seems to be what we call modulational instability p8[ , 
not rogue photodissociation. 

In sum, we have proposed and numerically solved a 
minimal but realistic model for the photodissociation of 
molecules to noncondensate atoms that invariably ac- 
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companies coherent photoassociation of a Bose-Einstein 
condensate. The results support the notion that, as the 
atoms need time to collide, there is a maximum achiev- 
able rate to photoassociation of a condensate. What is 
dimensionally obvious but might be physically surpris- 
ing is that the maximum rate in a noninteracting, zero- 
temperature condensate should be of the order of the 
only frequency parameter that can be constructed out of 
the density of the condensate. Apart from the fact that 
rogue photodissociation sets experimental constraints on 
coherent photoassociation, experiments would also be in- 
teresting because they bear on the wave versus particle 
nature of atoms, and on the limits of validity of scattering 
theory. 

This work is supported by NSF Grants PHY-9801888 
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carried out in part at the Aspen Center for Physics. The 
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photoassociation rate was first made to the authors by 
Wolfgang Ketterle. 
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FIG. 1. Atomic (Pa, dotted line) and molecular (Pm, 
dashed line) probabilities, as well as the total probability that 
an atom is part of either condensate (Pa + Pm, solid line) as 
a function of time. In panel (a) there is no rogue photodis- 
sociation, 7 = 0, while in panel (b) we set 7 = l/(27r). The 
(renormalized) detuning is fixed at 5 — £1. 
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FIG. 2. Maximum rate of photoassociation R for a fixed 
interaction time r scaled to the density dependent frequency 
parameter u p as a function of 1/(tuj p ). The lines are for 
various Rabi frequencies f2, as indicated in the legend. Phys- 
ically, in order to follow a line of fixed Rabi frequency toward 
increasing 1/(tu> p ), one must at the same time increase the 
laser intensity and decrease the condensate density, so that 
rogue photodissociation becomes more prominent. 
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FIG. 3. Time dependences of the three probabilities Pm 
(dotted line), Pa (dashed line), and total condensate proba- 
bility Pa + Pm (solid line) as a function of time. The fixed 
parameters are Q = 4/r, S — —9.27 0,, 7 — 1.27. A true 
photoassociation rate does not exist for these parameters. 
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